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We consider the Hubbard model and its extensions on bipartite lattices. We define a dynamical 
^^ ' group based on the 77-pairing operators introduced by C. N. Yang, and define coherent pairing states, 

which are combinations of eigenfunctions of 7)-operators. These states permit exact calculations of 
numerous physical properties of the system, including energy, various fluctuations and correlation 
Tx I functions, including pairing ODLRO to all orders. This approach is complementary to BCS, in that 

these are superconducting coherent states associated with the exact model, although they are not 
eigenstates of the Hamiltonian. 
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The Hubbard model plays a special role in condensed matter physics. It allows one, in appropriate limits, to model 
I . the electronic properties of systems ranging from insulators to superconductors. It is generally believed that high-Tc 
' ^ ' superconductivity may be described by some form of the Hubbard model. Although the model cannot be solved 
except in one dimension, some insight into its properties in general dimensions can be obtained through the so called 
?7-pairing mechanism introduced by C. N. Yang |1|. This mechanism allows one to construct a subset of the exact 
spectrum of the model. The eigenfunctions obtained through jy-pairing possess the property of off-diagonal long-range 
order (ODLRO) and thus are superconducting. In this note we introduce a new family of wave functions which are 
combinations of ?7-pairing eigenfunctions. The 77-pairing procedure has been applied to a number of strongly-correlated 
fermion systems B]- M. Our wavefunctions are coherent pairing states (CPS) of the dynamical group of the Hubbard 
rs^ model. Although not eigenfunctions of the Hamiltonian, they permit exact calculations of numerous physical proper- 

^Nj ties of the Hubbard model, including the energy, arbitrary moments of the Hamiltonian, fluctuations and correlation 

OQ functions, including ODLRO which is shown to be non-vanishing. The CPS are mathematically related to the varia- 

^~~^ \ tional wave functions used in a mean field treatment of BCS type H • 

i^ I ry-pairing For the Hubbard model we adopt the definition and notation of Yang m. Let ai and bp be real-space 

j^ ' creation operators for spin- up and spin-down electrons respectively, i.e. ci= ai , ci = bi with ai and bi satisfy- 

Ch ] ing the usual fermion anti-commutation relations. 

I ■ Consider a 3D Hubbard model on.a,LxLxL = M cube (L even) with periodic boundary conditions. The 

■ O ' Hamiltonian is given by 

O. H = n+Ti + V (1) 



To = AeY^{aiaj. + bib^) (2) 

k 

Ti = — B y^icos kx + cos ky + coskz){aia~ + bibj-) (3) 

fe 

V^2Wj24arbpr, (4) 
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where e > 0, ai is the Fourier transform of at , 2W is the on-site Hubbard interaction of arbitrary sign, and A and 
B arbitrary constants. We introduce the ry-operators which create (annihilate) a fermion pair with momentum tt: 



Y^ e^^aphp 



Y. H^^-k^ 



(5) 



n" -Y.^''4^ -T.K-A- 



(6) 



It has been shown ( H], @) that the operator 7^+ satisfies 

[H,rj+] =Erj+, 



(7) 



with E — 2Ae + 2W. Equation ( 0) is typical of a Spectrum Generating Algebra ( 10, |Tj]) and implies that for any 
power-expandable f(r]~^) 



[HJ{rj+)] =E^+n^^ 



(8) 



Note that E does not depend on B. The relation (^) for the Hubbard model was derived some time ago Q but its 
consequences were only fully exploited by Yang p]. The operators ?/ satisfy the angular momentum commutation 
relations of SU(2): 



[77+, 77] = 277^ 



(9) 



1 Y^/ (O) , (b) ,^ 

= 2 2^^4+4 -1) 

r 



(10) 



where the local occupation number np is equal to rip ^ rip — ata^ -f btbp. We also observe that from Eq.(p|) that 
the following relation holds 



r] r]^ 



M VM 



= l-d, 



(11) 



where d is the electronic density, d ~ M^^{^pnp). Equation (^l|) indicates that for small electron density the 
operators -p= are approximately bosons [|2|. The operators ry also satisfy the relations (77+) = (r/) =0, 

reflecting, according to the Pauli-principle, the impossibility of occupying a given site f by more than one pair 
(a^6^). For given M and using (7) one can produce M exact, normalized eigenstates of H by applying succesive 
powers of ?/+ on the vacuum state |i;ac). So 



|*Ar)=/3(iV,A'f)(r;+)"|mc), N = 1,...,M 



(12) 



is a simultaneous eigenstate of H and of the operator N2 counting the number of doubly occupied sites, N2 

V^ (a) (6) 

z ^r r r ' 



where (3{N, M) is a normalization factor equal to ||l[ 



m,M) 



{M -N)l 
M\N\ 



(13) 



(14) 



Evidently, {'^n\{ii^Y\'^n) = Sr.o- Note that [H,N2] 7^ 0. We observe that the \'^n) depend neither on the value nor 
on the sign of W. In general, [77, H] j^ except for the half- filled band IQ. 



Dynamical Group for H With this in mind, we embed the Hamihonian H together with ?/= {?/, ^7 , ?yz}, in a- larger, 
dynamical, group. Define a new operator Jq by 

H , , 

Jo = ^-Vz- (15) 

Using (M) and its hcrmitian conjugate we find that 

[Jo,rf\=0. (16) 

We conclude that the smallest group containing H is {Jo,fj}, where Jq is the center of the group but not the unit 
operator. The dynamical group of our Hubbard model is thus U(2). This would appear to be the first instance of a 
dynamical group for an exact interacting many-body system. The relation (|l^) is essential for the calculation of any 
expectation values of H. 

Coherent Pairing States We introduce a normalized spin coherent state by 

= il + \tA')-^e'^^\0) (17) 

0) = -pi-(?7^)*^|vac). We refer to |/i > as a coherent pairing state. This step 
q|, which is however not related to any Hamiltonian with a local potential 



where the state |0) is the filled pair state 

is reminiscent of the BCS wave function 

energy. In contrast, our states arise out of the exact relations Eqn.(|7|) and Eqn.(16). In the limit M -^ oo, |/i) 

becomes an eigenstate of -p= and apart from normalization is a harmonic oscillator coherent state |14 . The state 

l/i) is not an eigenstate of H. In contrast with J^at) it involves components with different numbers of particles (pairs) 

and thus gives rise to non-zero values of {fj,\r]^\ii). Further, using (|7|), (P) and (|l6|) we may calculate (/i|iJ^|/^) for any 

p = 1,2, 3, ... in terms of {p,\ri^\p). We first calculate (/i|i?|/i) by purely algebraic means: 

H\^i)^N-i{[H,e^^]+e'^^H}^irj+)''\vac) 

= A/'-3 {^Eriiie^"'^ + e^^'^ME} —^{r]+f'^\vac) 
= {-liEi] + ME)\ii). 

where J\f{n) = (1 -|- |/ip)*^. The required expectation value {fi\H\fi) becomes AIE — ^iE{^i\ri\^i), which, using the 
results of [Q (Formula 4.2), leads to 

, , , > ME , , 

The formula (jlj) indicates that the energy of the state \\i) (which involves different numbers of pairs) is equal to 
the energy of the fully filled state |0) (ME) reduced by the factor (1 -I- |/ip)~^ < 1. The physical meaning of the 
parameter p, is obtained from the average number of pairs in a state \pi), (/i|iV2|/i). Since |/i) does not depend on the 
Hamiltonian's parameters, 

, ^.r ^ ^ ^ , ^^il ^ ^ ^ M dE 

(.IW = -(M|^|M)^^^^^^. (19) 

We conclude that 

H' = ^-i, (20) 

where n-i — ^-^ is the average density of pairs in the state |/j,). We may extend the set of states for which exact 
analysis is available by introducing r-depleted states, defined by (normalized) 

\y.-r)^Mr-'^rf\pi). (21) 



These are analogues of the displaced number states of quantum optics [|5| . These states give rise to a more interesting 
energy spectrum than the equidistant Yang case, with the gap between neighbouring depleted states |/i; r) and |/z; r— 1) 
being given by 

^ n |2^ _ {^l\(7l+r-'Hln'-^\^i) {^,\{rJ+YH,^^■\^^) 

For /i = we evidently have Yang's functions for which all the gaps are strictly equal to E. For fJ. j^ we expect a 
structure in A^d/ip). In fact all the quantities in ( P4 ) can be calculated using only A/'(|/x| ) and (/i|i/|/i) ([iq). For a 
general operator Q we can calculate (m|(?7^)'^Q'7'^|m) through the relation 

(Ml(r?+)''Qr;'-|M) = (1 + ■/^^"'''^(^^ [(1 + iMH^'iMlQI/^)] , (23) 

which indicates that for p = |//p the generating function for the matrix elements of Q between the depleted states is 
proportional to (1 + p)*^(/i|Q|u}, which for Q = 1 and Q = H furnishes all the input for Eq. (p2[). 

The detailed analysis of Eq.(E2h confirms a very interesting structure of A^ as a function of p, M and r. The precise 
description will be presented elsewhere but we note here that the gaps as a function of p go through a maximum for 
r « ^ which in turn disappears for r > ^. This confirms that the half- filling point {N = ■^) plays a special role 
for the Hubbard model. Equation ES may be used to obtain the following simple result for the energy dispersion in a 
coherent pairing state: 

(^^)^ - ^ (24) 



{p\H\p)^ M 

where p — |/zp. This indicates that the energy fluctuations are normal in the thermodynamic sense, as in the grand 
canonical ensemble. Similarly, in the first depleted state \p; 1) 

(Ag)i' ^ p(2 + 2Mp - Mp^ + APp^) 
(F)i2 {M -l){l -p + Mpf ■ ^ ^ 

Note that the dispersion in the first depleted state Equation(p5|) is always greater than that in an SCS, Equation (p4|). 
Analogous, if more complex, results hold for higher depleted states. 

Time evolution of the Coherent Pairing States Since the coherent pairing states are not eigenstates of the Hubbard 
hamiltonian, they possess a non-trivial time dependence. This time evolution is easily calculable via the time- 
dependent Schroedinger equation due to the simple algebraic structure of the model. For the case of a conventional 
coherent state satisfying a\z >— z\z > evolving under the action of a hamiltonian H = ijj{a^a + ^), the evolution is 
simply expressed by the propagator {h = 1) 

I < z{0)\z{t) > I = exp{\z\^[cosujt - 1]). (26) 

In the case of the coherent pairing state Eq.(ll7|), the analogous result, with Eq(|7|) is 

l<MO)Kt)>H|^^f|L_^|. (27) 



In the limit Af — > cx), /i — > — ^ (corresponding to the group contraction rj —^ yMa, compare Eq.([ll|) ) we recover 



the conventional (bosonic) case (|2f 

Off-Diagonal Long-Range Order (ODLRO) The presence of Off-Diagonal Long Range Order (ODLRO) iQ is de- 
tected by the nonvanishing of correlators such as < a^fo^forOr > as |r — s| — > 00. Yang has shown that his states display 
ODLRO which, in the thermodynamic limit, is proportional to 712(1 — 71.2), where n2 = N/M is the pair density. We 
may similarly show that our SCS states |/Lt), /x 7^ 0, exhibit ODLRO, also proportional to 7X2(1 — ^2) where the average 
pair density 7i2 =< N > /M. Additionally, the states \p;r) exhibit ODLRO and all the results reduce to those of 
Yang ior p — [O. Thus the states \p) , \p; r) are superconducting for all p and r. It is worth noting that although 
{ip N\'n\''P n) = , which makes 77 unsuitable for defining an order parameter in the usual sense, {p.\ri\p) 7^ as in the 
analogous BCS case. 



Relation to Mean Field Theory We may now write a mean- field version of the Hubbard Hamiltonian H — J2k-^i 



k 



Hk = Ekiak^ak + bk^k) + 2W{^k*m + ^kVk^) (28) 

with rjq ~ ^j, akhq-k {q ~ tt) and effective energies Ek which include the Tq and Ti terms of Equations(g) and (|^). 
The spectrum-generating algebra for the Hamiltonian Eqn(28) is Ufc(2). The non-diagonal terms are the number- non- 
conserving an alog ues of a spin-density wave system (see, for example, |L8| ) and were already observed in a multi-phase 
SU{8) model |1^, where they were called "anomalous" terms, their relation to the Hubbard model at that time not 
having been appreciated. Thus in the mean-field approximation the dynamical group is 'S)kUk(2). The associated 
group parameters (Bogoliubov transformation angles) are /ik, vector analogues of the /i parameter in the U(2) which 
is what remains of the dynamical symmetry in the exact model. 

Extensions beyond the Hubbard Model It is possible to demonstrate that the relation (^ (with modified E) can 



be fulfilled by Hamiltonians other than (1). Yang already mentioned the possibility |l[| of non-local interactions 
satisfying (M) with appropriately modified rj. It is interesting to observe that at least one case of a truly non-local 
interaction satisfies (0). It concerns an extension of the pair- hopping model [^ of the form ||l^] 

Hph = To + Ti + F^ 77; • 77; - - ^ n, (29) 

(rs) r 

which satisfies [Hph — (Tq + Ti), 77+] = Vr]^. Still other extensions are possible for which the coherent pairing state 
1^) is a useful tool pTf. 
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